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ABSTRACT. In this paper, we attempt to describe the spread of COVID-19 based on an 
extensions of the well-known SEIR models. We bring in a new representation to appraise 
and manage the spread of infectious disease COVID-19 through SEIHR pandemic model, 
which is based on the supposition that the infected individuals are isolated and the necessary 
treatments are arranged so that they cannot taint the other residents in the community. 
Dynamics of the SEIHR model is presented by basic reproduction number Ro and the 
comprehensive stability analysis. The reproduction number was also determined for the 
epidemiological model and found to be consistent with the local stability condition for the 
disease-free equilibrium. We establish the positivity and boundedness of solutions, local and 
global stability analysis of equilibria to examine its epidemiological relevance. To validate 
the model and estimating the important model parameters and prediction about the disease, 
we consider the real cases of Ethiopia from 4th March to 20th August 2020. 


1. INTRODUCTION 


In December 2019, corona virus first reported in Wuhan, China, is an infectious disease 
caused by a newly discovered corona virus. The virus is now known as the severe acute 
respiratory syndrome corona virus 2 (SARS-CoV-2). The disease it causes is called corona 
virus disease 2019 (COVID-19) [1]. According to WHO Novel Corona virus disease situa- 
tion report 21 January 2020, the first case outside of China was reported in Thailand on 
13 January 2020, Japan on 15 January 2020 and Republic of Korea on 20 January 2020 
[2]. Currently the COVID-19 outbreak has now spread over the world. According to WHO 
the symptoms of COVID-19 may appear two to 14 days after exposure. This time after 
exposure and before having symptoms is called the incubation period i.e., during this period 
the infected humans are not infectious. However, the patients infected by COVID-19 will 
have symptoms including fever, cough and difficulty breathing, and the severe patients may 
have renal failure. This list is not all inclusive [3]. Current evidence suggests that COVID-19 
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spreads between people through direct, indirect (through contaminated objects or surfaces), 
or close contact with infected people via mouth and nose secretions. These include saliva, 
respiratory secretions or secretion droplets. These are released from the mouth or nose when 
an infected person coughs, sneezes, speaks or sings [1]. According to WHO, to prevent in- 
fection and to slow transmission of COVID-19, do the following: wash your hands regularly 
with soap and water, or clean them with alcohol-based hand rub, maintain at least 1 me- 
tre distance between you and people coughing or sneezing, avoid touching your face, cover 
your mouth and nose when coughing or sneezing, stay home if you feel unwell, refrain from 
smoking and other activities that weaken the lungs, practice physical distancing by avoiding 
unnecessary travel and staying away from large groups of people. According to WHO 2020 
report, after the outbreak day to present September 18, 2020 globally, there were about 
30,055,710 confirmed COVID-19 cases with 943,433 deaths. In Africa also the virus spread 
to 52 countries, but two countries (Lesotho and Comoros) are still registered as COVID-19 
free. Through out Africa up to August 26, 2020 confirmed cases are 22,313, number of deaths 
are 1,124 and recoveries are 5,492 [4]. In Ethiopia the first case was reported on March 13, 
2020, and up to September 18, 2020 the confirmed cases reached to 67,515, deaths are 1,072 
and total recoveries are 27,638. In epidemiology, mathematical modeling has a great role 
in understanding the dynamics of infectious diseases in a population and widely used to 
predict the results of an epidemic successfully. Most commonly used epidemic models are 
SIS, SIR and SEIR models [4]. The SIR model is a very well established model and used 
widely for various epidemics [5]. Several Mathematical modeling have been proposed to 
study the dynamics of epidemic disease. Mathematical models that analyze the spread of 
COVID-19 have began to appear in published papers and online resources [6, 7, 8]. A num- 
ber of compartmental models have already been proposed and analyzed for the COVID-19 
outbreak in different countries {9]-[19]. For instance, [9] included compartments correspond- 
ing to suspected cases, which consists of the individuals that show similar symptoms but are 
not confirmed cases, and indirectly infected individuals. The modified SEIR model which 
included asymptomatic and treatment compartments for occurrences in Wuhan, China, the 
city where the outbreak started, and outside of Wuhan haven been used [10]. [12] included 
a separate compartment for quarantined individuals in the SIR model to account for the 
containment measures applied by the public for the virus. The stability of epidemic models 
has been studied in many papers [5]-[22]. Many authors paid attention to local stability of 
equilibria. Recently, the study of epidemic models mainly concerns global asymptotic stabil- 
ity. The most successful approaches to the problem are the direct Lyapunov method [23] and 
the geometric method [24]. In this paper, we consider a modefied SEIHR model depicting 
the transmission of corona virus with disease resistance in human. The model is given by 
a system of four differential equations depending on parameters. By using the method of 
next generation matrix [25], we found a threshold Ro called basic reproduction number. In 
general, when Ro < 1, the disease dies out and when Ro > 1, the disease persists in the 
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population. If we suppose that the endemic equilibrium also exists for Ro < 1, although it 
is not true, then the bifurcation occurring in the model can be explained as a trans critical 
bifurcation. Several various methods are used to determine the stability of equilibria. We 
concentrate our study on the globally stability of equilibria. This is obtained by Lyapunov 
functional approach and geometric approach 


2. MODEL FORMULATION 


In this section, we develop mathematical model for the transmission of COVID-19 which 
spreads in a populations. Based on the development and epidemiological characteristics of 
COVID-19 infection, we extend the classical SEIR model to describe the transmission of 
COVID-19 in Ethiopia. In this model the entire population is considered into four disjoint 
compartments. Susceptible individuals (S) are those who are not infected by the disease 
but there is a possibility to be infectious. Exposed individuals (E) are individuals who are 
in the incubation period after being infected by the disease, and have no visible clinical 
signs. This individuals could infect other people with a higher probability than people in the 
infectious compartments. After the incubation period, the person passes to the Infectious 
compartment; Infected Individuals (I) are individuals who developed the symptom of the 
disease. As far as the diagnosed individual are concerned, we assume that infected individuals 
will be immediately sent to designated hospitals for isolation and treatment, these group of 
the population will be converted to hospitalized H. Finally, Recovered individuals (R) are 
those individuals who recovered (through treatment or natural recovery) from the disease. 

The recruitment rate denoted by 7 is either through immigration or birth. The susceptible 
individuals got corona virus disease by contact rate of @ either from exposed individuals 
with probability of 1 — 7 or from infected individuals with probability of r and move to 
the exposed compartment. The exposed individuals become infectious and join the infected 
compartment with the proportion a. ø is the average rate at which symptomatic individuals 
become hospitalized and e€ is the recovery rate without being hospitalized. y is the recovery 
rate of the isolated infected by treatment from the disease or die due to the corona virus 
disease at a rate of 6. The recovered individuals become again susceptible to the disease 
with a rate of 7. 

Based on the above assumptions and the actual isolation strategy, the spread of COVID-19 
in the populations is depicted in Figure 1. 
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FIGURE 1. Flow diagram of the compartmental model of COVID-19 


Using the above depiction, we formulate the corresponding dynamical model as follows. 


dS 

ap = TA + (1-1) EIS +R — pS 

dE 

ET = p|(TI +(1-T)E|S — (a + u)E 
(2.1) d =a — (o +e +p) , 

dH 

dR 

ap UR 


With the initial condition 


There are a total of 12 model parameters in system (2.1). 





Parameters | Description 





T Recruitment rate of individuals 

B Contact rate of individuals 

T Transmission rate of infection 

a Proportion of exposed individuals leaving the compart- 


ment 

recovery rate without being hospitalized 

Proposition of infected individuals become hospitalized 
Recovery rate of hospitalized patients 

Proportion of recovered individuals to be susceptible 
Natural death rate 


Death rate due to coronavirus 
TABLE 1. Description of the model parameters 
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3. ANALYSIS OF THE MODEL 
In this section, we study the quantitative and qualitative analysis of Eqn. (2.1). 


3.1. Well-posedness. An important feature of an epidemiological relevant model is the 
positivity and boundedness of the solutions. Therefore, it is important to prove that all the 
variables are nonnegative for all time t > 0 which implies that any solution that has positive 
initial values will remain positive for all time t > 0. We begin by determining the biologically 
feasible set for the model (2.1). 





Theorem 3.1. The closed region Q = {(S,E,I,H,R) € RË : N(t) < at is positively 
invariant set for the system (2.1). 











Proof. For the model the total population is 
N = S(t) + E(t) + I(t) + H(t) + R(t) 





Then, differentiating N with respect to time t, we obtain: 
dN dS dE d? dH dk 





= | | | H = = ôH — uN 
dt d dd dt dd d g 
If there is no death due to the disease, we get 
dN 
el ae ae nail 
(3.1) ap 


Integrating the inequality 3.1, using initial condition, we obtain 
N(t) < N(0)e™ + = [1 — e=] 
u 


Evaluating it as t > oo, we get 
N(t) < 


Therefore 





Q= {(S,E,I,H,R) ERË : N(t) < 








} 


Which is the feasible solution set for the model (2.1) and all the solution set is bounded in 
it. 





RIN 














3.2. Asymptomatic Stability of Disease Free Equilibrium. When there is no disease 
in the population, I.e E = I = 0, the disease free equilibrium occur and to determine the 
equilibrium solutions, we set the right side of Eq. (2.1) equal to zero and obtain 


r—Bi(rTI+(1—T)E|S+nR-pS = 0 
Bl(7l+Q1-7TNE|S-(a+pe = 

(3.2) aE —(o+e+pu)I = 
ol-(y+é6+p)H = 

el+yH-(ntWR = 


D. = a ee 
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Therefore the disease free equilibrium point is obtained to be 
T 

(3.3) £9 = (=.0,0,0,0) 
u 


Basic reproduction number. The basic reproduction number is an important threshold con- 
dition in the analysis of an infectious disease. It determines whether the disease will die out 
or persist in the population as time increases. We calculate the basic reproduction number 
Ro of the system by applying the next generation matrix method used in [?]-[11]. The first 
step to get Ro is rewritten the model equations starting with newly infective classes: 


dE 

= Bll +1- TES — (a + p)E 
(3.4) Z =ak—(o+tet+pyl ; 

dH 


Then by the principle of next-generation matrix, the Jacobian matrices at DFE is given by 





Plo orn) at fi 0 0 
F= 0 0 0|, V= —& o+éet+p 0 ; 
0 0 0 0 —0 y+ô+u 
Then 
B-—r)r(o+e+u)+8TTra BTT 
ula+u)(o+e+u) ulo+e+n) 
FV“ = 0 0 o|, 
0 0 0 


Therefore, the basic reproduction number is given us 


P1- T)r(o +e + u) + Braa 


oe CO ome cep) 





which is is a threshold parameter that represents the average number of infection caused by 
one infectious individual when introduced in the susceptible population [19]. 


Local Stability of DFE. Local stability is attained when the eigenvalues of the Jacobian, 
A, are negative or have negative real parts. In other words, the solutions for A such that 
[Jæ — AZ| = 0 should be negative or have negative real parts if the solution is complex [27]. 


Theorem 3.2. The DFE £o of the model (2.1) is locally asymptotically stable if Ro < 1 and 
unstable if Ro > 1. 
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Proof. The Jacobian matrix at £ọ is as: 


B(l=7)r Brr 


TH a © p ” 

0 AD ay frr 0 
Ja=ļ| 0 q —0 —E€— pH 0 0 

0 0 a —y-oéd-p 0 

0 0 E 7 Say 


and the characteristic polynomial is |J-, — AZ| = 0. Solving this polynomial, we get: 











(3.6) (=u — A) (=7 — 6 — u — A) (=n — u — A) (X + AA+ B) =0 
where 
ES 
r SDRE U sip Sey Bd 
u 
B —rT)r BTTQ 
B EE LUNT (a+u) (o +e +u) 
H H 
Thus, the eigen values of the first expression of eqn. (3.6) are given by 
Al = —u<0 
à2 = —(ņ+p4)<0 
AB = -(y+d+p) <0 





The last expression of eqn. (3.6), that is 
(3.7) + AA+B=0 


We applied Rourth-Hurwitz criteria and by the principle eqn. (3.6) has strictly negative real 
root iff A> 0,B > 0 and AB > 0. Clearly, we see that A > 0 because it is the sum of 


positive parameters. 


fi a otetu BTT 
atu  pla+u)(o+e+u) 


B = (a+ p)(0+e+ u)(1-— Ro) > 0 





I= Ro>0 











Therefore, £o is locally asymptotically stable. 





Global Stability of DFE. 


Theorem 3.3. The DFE £o of the model (2.1) is globally asymptotically stable if Ro < 1 
and unstable if Ro > 1. 


Proof. To prove the global stability of the equilibrium point we construct the Lyapunov 


function as 


(3.8) L= GE + Col 
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Differentiating equation (3.13) with respect to t gives 
dL dE dl 

















(3.9) dt = Cy dt T C2 di 
Substituting @ and # from the model (2.1), we obtain: 

dL 

ne (BITI + 1 -T)EIS — (a + HJE) + (ak -—(o+e+p)l) 
(3.10) = apl- rT)ES -— cla + u)E + caE + caprI — clo +e +u 

Take cy = SETO then we get: 
dL 1— — 
(3.11) ae GbTIS — (-4 oe Cida) (otet+p)l 
(3.12) = Abate sete ee eee eet) ii 
a 


Taking cp = 1, and substituting Ro we get 
dL fi ad, 
dt 7 ua 
2 (£ uo +e #2) (e e aes -1) j 
pla + p)(o +e + u) 


= (£ ula +e P) (Roi) 
































for S < So < 7 and dL < 0 for Ro < 1 and trajectory of the model (2.1) on which aL = 0 if 
and only if J = 0. This implies that the only a < 0 is €9. Therefore by Lasalle’s invariance 


principle, €o is globally asymptotically stable in Q. 














The Endemic Equilibrium point. The endemic equilibrium point shows that the disease will 
persist in the system in the steady state. Here solve the equations (3.2) to obtain S, E, I, H 
and R. But for easy identification S, Æ, I, H, R are represented by (S*, E*, I*, H*, R* ) at 


the steady state of the endemic respectively. 























g — Mateo +e tp) 
Blrat+ (1—-T)(o+e+p)| 
pe Leten whe ts H u)’ (y +ô +H) (Ro—1) 
(3.13) r eu in E (Ro —1) l 
yt = Orla wha ie be +h) (Ro —1) 





pr — SL + win + wo +e + ulel +ô + u) + 72] 
BK 





(Ro — 1) 
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where 





k = |ra+(1—-rT)(ot+et+p|[(atw(nt+p(otetp)(yt+ot+p)—na((yt+6+p) +70) 


Theorem 3.4. The endemic equilibrium point eğ of model (2.1) is locally asymptotically 
stable in Q if Ro > 1. 


Proof. Recall the Jacobian of the system (2.1) at any equilibrium point (S, E, I, H, R) and 


we have 


—p|rI + (1 -7)E] — u —B(1—7)S —BrS 0 n 
BITI + (1-r7)E] Bl-Tr)S-a-p prs 0 0 
J= 0 Qq —0 —E— pH 0 0 
0 0 og —y—-ô-—u 0 
0 0 E 7 =p 


At the endemic equilibrium point ej, calculating the Jacobian matrix J and then solving 
det(J — AI), the characteristic equation is 


(3.14) (—y=6 = p—A)(—9 — u — A)(X? + aà? + DA+c) 
where 
a = B(1-7)S*+8(TI+(1-7)E)-34-a-0-e 
b = [=g -r7)S* +2u+a+o +e] l8(TI+(1-rT)E)-— u] 
— (6 — T)S*B(TI + (1 -—7)E)+ (a+ u+ e+ p)| 
c = B7rI+(1-7)F)[(6+¢e+4+4)(80. —7)S* + (a+ u))-— Bras*] 
—uļla +u) (ô +e + p) 








Since roots of the characteristic equations corresponding to J(eġ) are 


à = —(y+6+p) <0 
Ag = (n+p) <0 





and other three satisfies the following cubic equation 
N+ ad? + bA+c=0 


applied Routh-Hurwitz criteria. 














4. SENSITIVITY ANALYSIS 


We carried out sensitivity analysis, on the basic parameters, to identify their effect to 
the transmition of the disease. To go through sensitivity analysis, we used the normalized 
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sensitivity index definition as defined in [29]. The Normalized forward sensitivity index of a 
variable, Ro, that depends differentiably on a parameter, u, is defined as: 

= OR» u 


Ar 
k ðu ERa 


for u represents all the basic parameters. Here we have 


BA-—T)r(o+e+p)+ brra 





















































R p= 
; (a+ u)(o +e + p) 
For the sensitivity index of Ro to the parameters: 
ORo | B 
ABR se sisu 
3 T x R > 
ORo T 
Ni a eg ee a | 
a r x R, >0 
Ro ORo T Br(o+tet+yp) 
APL = = x — =l- 
ðr Ro ato + E+ wo 
Ro . Ro a Brulote+nlr-(-7)(o+e+ n) 
Aa S x = <0 
Ba * Ro Mat wo+e+ uP 
ARo = ORo g&n Brraula + u) 22 
de “Ry [nla +uy(o +e +p) 
ARo ORo TOO Brraula + u) 22) 
ðo Ry ulate be +p 
2 
ARo ORo U noa u)(o H) — 9 
On Ro wat p)(o H) 











The sensitivity indices of the basic reproductive number with respect to main parameters 
are 7,9 and T. Which show that they have great impact on expanding the disease in the 
community if their values are increasing. Also those parameters in which their sensitivity 
indices are negative a,€, and u have an effect of minimizing the burden of the disease in 
the community as their values increase. Therefore, research advice for steack holders to work 


on decreasing the positive indeces and increasing negative indices parameters 


5. NUMERICAL SIMULATIONS 


In this section, we examine the COVID 19 model and studied the effects of combined 
strategies on controlling the transmission of the disease. We start to solve the state equations 
with a guess for the controls over the simulated time using the fourth order Runge-Kutta 
scheme. Because of the transversality conditions, the adjoint equations were solved by a 
backward fourth order Runge-Kutta scheme using the current iterations’ solutions of the 
state equation. This process is repeated and iterations are stopped if the values of the 
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x10’ Predicted Spread of COVID-19 in Ethiopia 
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unknowns at the previous iteration are very close to the ones at the present iteration. For 
numerical simulation purpose, we have used empiricaldata of COVID-19 cases in Ethiopia. 


Moreover up to the April 23, 2020 total recoveries were 21 and registered deaths were 3. 
But, after march 20, 2020, COVID-19 in Ethiopia. 


6. CONCLUSION 


In this current paper we have formulated and studied an epidemic model of COVID- 
19 virus which is transferred from human to human. Therefore, prediction about infected 
individual is very much important for health concern arrangement of the citizens. It is also 
important to control spread rate of the COVID-19 virus with restricted supply. In this 
study, we have formulated a SEIHRS epidemic model for pandemic COVID-19. Considering 
all major parameters of the progression of the disease, several analytic results established. 
Theoretically it is proven that the dynamics depends on the basic reproduction number 
to examine the stability of the system. All the properties necessary for epidemiological 
relevance have also been proved. We have estimated the parametric values for Ethiopian 
using the existing real discrete data. The acceptable agreement between data analysis and 
numerical solutions are established. Finally the inclusion is that to estimate the situation of 
global pandemic COVID-19, mathematical modelling is an efficient method if the parameters 
can be estimated properly. We propose a novel dynamic system with time delay to predict 
the trend of outbreak for the 2019-nCoV. In this model, the instantaneous increment of 
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cumulative diagnosed people depends on the history of cumulative infected people, by which 


the latent period can be taken into consideration. The numerical simulation and parameter 


identification were carried out to verify the effectiveness and accuracy of the novel dynamic 


system. In addition, this model can well approximate the true data in this event, which may 


further predict the tendency of event. 
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